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ANALYSIS  OF  EXPERIMENTAL  THERMAL  VALIJES  AND 


EQUATION  OP  STATE  FOR  HYDROOEN 
V.  A.  Rabinovich 


Detailed  analysis  Is  made  of  the  p  v  —  T  data 
of  various  experimenters  and  this  Is  used  to  derive 
a  grid  of  reference  values  for  thermal  quantities. 

An  equation  of  state  Is  set  up  for  hydrogen  and  the 
feasibility  of  extrapolating  It  Into  the  hlgh»tempera- 
ture  region  while  allowing  for  the  dissociation  of 
hydrogen  molecules  Is  demonstrated. 


The  majority  of  experimental  data  on  the  thermal  properties  of 
substances  are  obtained  In  the  best  case  from  the  curve  of  satxxra- 
tlon  IQ)  to  a  teniperature  on  the  order  of  600®C.  The  experimental 
caloric  values  In  the  literature  are  relatively  scanty  and  In  the 
majority  of  cases  are  completely  lacking. 

The  problem  of  the  present  work  was  to  derive  an  accurate  equa¬ 
tion  of  state  for  hydrogen,  correct  In  the  gaseous  phase  In  a  broad 
Interval  of  change  In  temperature  and  pressure. 

Analysis  and'  Concordance  of  Experimental  and  Thermal  Data 

It  Is  necessary  to  create  a  network  of  grid  reference  data  on 
existing  thermal  values  In  order  to  set  xq;)  an  equation  of  state  and 
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to  disclose  the  patterns  and  peculiarities  of  the  thermodynamic  behav¬ 
ior  of  a  substance. 

The  Initial  values  were  graphically  conformed  In  dimensionless 
coordinates  along  the  Isotherms  a-iT  >  f((D)  and  the  Isochores  o-r  • 

■  f (t)  ,  where  p  -  pv/llTj^,  a>  ■  Vj^/v,  and  t  - 

As  a  result  of  this  concordance  of  a  large  number  of  experimental 
p  -.i-  V  — >  T  data  [1-9]  we  were  successful  In  analyzing  and  estimating 
the  accuracy  of  the  Initial  values. 

In  the  analysis  It  was  discovered  that  the  data  In  Bartlett  [1] 
and  Wlebe  and  daddy  [2]  differ  from  each  other  by  vq;)  to  but  their 
average  values  link  up  well  with  the  experiment  In  Michels  et  al.  ['}] 
idilch  Is  distinguished  by  a  high  accwacy  noted  by  the  authors  them^ 
selves  and  confirmed  on  constructing  Isothermal  and  Isochorlc  cross- 
sections.  The  data  In  Holborn  and  Otto  [4]  are  slightly  Inferior  In 
accuracy  to  the  thermal  values  In  Michels  et  al.  [3]«  tut  within 
limits  of  Q»2^  they  may  be  smoothly  linked  with  the  latter.  The 
ejqperlmental  valuea  In  WltkowslcL  [3]  and  Onnes  and  Braak  [7l  are 
excessive  and  with  the  growth  of  o>  the  devlatl<»i8  frcm  Michels  et  al. 

[3]  and  Holborn  and  Otto  [4]  reach  0.^.  The  data  In  Onnes  and 
Fenlng  [6]  also  completely  satisfactorily  conform  In  the  low-tenpera- 
ture  regions;  at  temperatures  near  the  critical,  however,  noticeable 
discrepancies,  reaching  2$  close  to  the  critical  point,  are  observed. 

The  siq|>ercrltlcal  point  was  therefore  the  object  of  the  especially 
careful  processing  of  the  ejqperlaental  data. 

The  graphic  concordance  of  the  experimental  p  —  v  —  T  data 
enabled  tis  to  create  a  grid  of  thermal  reference  values;  with  this  as 
a  basis  we  were  successful  In  disclosing  several  patterns  In  the  thnmo- 
dynamlc  behavior  of  hydrogen  cosqprlsed  In  the  following. 
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1.  All  the  Isochotres  In  the  experimental  Interval  of  tempera¬ 
ture  chaige  0.6  to  20  with  o)  <  1  have  an  upward-directed  convexity j 

l.e,,  <  0.  The  curvature  of  the  laochoree  Increases 

o  ■  0.4  to  0.6j  but  then  drops  anew  as  the  density  approaches  the 
critical. 

2.  In  the  range  a>  >  1  to  1.8  all  the  Isochores  have  a  double 
curvcture  and  the  point  of  Inflection  Is  the  region  of  change  r  «  1.7 
to  2.1.  When  eu  >  1.8  and  t  >  2.1  the  Inequality  (d*o/&T®)^  <  0  again 
becomes  true  for  all  Isochores. 

3^.  When  we  examined  the  Isochorlc  cross-sections  we  discovered 
no  stralght-llne  Isochores.  The  critical  Isochcn'e  Is  a  curved  llne« 
the  curve  of  which  decreases  as  It  approaches  the  critical  point. 

We  must  note  that  all  these  conclusions  are  In  complete  con¬ 
cordance  with  those  of  A.  Ye.  Sheyndlln  [10]  derived  from  the  presence 
of  maxlmums  of  c^  on  the  Isobars  and  the  anonalous  transitions  con¬ 
nected  with  them  In  pure  substances. 

Equation  of  State  for  Hydrogen 

The  best  known  equations  of  state  for  practicable  gases  (Beatty- 
Brldgeman,  Benedlct-Webb-Rubln«  Vtikalovlch-Novlkov,  and  the  one  set 
up  by  Ya.  Z.  Kazavchlnskly <  s  method  [11])  define  xn^esBure  £  as  a  func¬ 
tion  of  specific  volume  v  and  temperature.  T.  In  the  general  case  all 
these  equations  may  be  written  In  the  form 

P^fx(v)  +  h(f>)T  +  F(r,v).  (1) 

This  form  of  the  equation,  consisting  of  the  linear  portion 
ft(v)  fa(v)  T  and  of  a  function  composed  of  two  variables  F(T,  v) 
caused  by  the  distortion  of  the  Isochores,  proceeds  from  the  general 
thermodynamic  relatl<»i8hlps.  As  Investigations  wi  a  whole  series  of 
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•ubstanoes  hawshowni  function  F(T«  v)  is  small  In  comparison  with 
the  linear  portion  of  Eq.  (1)  and  monotonlcally  decreases  with  tempera¬ 
ture  Increase. 

We  selected  Kazavchlnskly's  equation  [11] «  which  derived  Its 
theoretical  basis  from  Kessel'man  [12]  and  In  the  general  case  Is 
written  In  dimensionless  parameters  as  follows: 


f 


to  describe  the  thermodynamic  properties  of  hydrogen.  Differentiating 
with  respect  to  r  we  obtain 

+ 7 (3) 

Rabinovich  [14]  has  shown  that  If  we  use  reliable  values  of  the 
curvilinear  portion  we  may  set  iq)  an  equation  of  state  which  very 
accurately  describes  not  only  the  thermal  values «  but  also  the  caloric. 

A  reference  grid  of  thermal  values,  from  which  ten  basic  Isochores 
were  selected,  was  used  for  separating  the  curvilinear  part  of  equation 
of  state  (2) . 

The  entire  curvilinear  portion  In  the  experimental  range  of 
hydrogen  parameter  change  was  obtained  In  tabular  form  as  a  result 
of  twofold  graphic  differentiation  with  Intermediate  and  final  flat¬ 
tening  of  the  lsoch(»*es. 

As  was  ascertained,  all  the  Isotherms  of  (d^o/dr^)  > 

(Fig.  1)  above  the  Boyle  temperature  T  >  comply  with  a  simple 
proportional  law,  l.e.,  the  ration  on  any  two  Isotherms  remains 

unchanged  vAien  d^  Is  the  sane. 

The  calculatl(ms  c<»iflrmed  that  the  entire  curvilinear  portion  of 
the  equatl(»i  of  state  for  temperatures  above  that  of  Boyle  may  be 
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described  by  an  analytical  expression  of  the  form 


at* 

The  method  of  determining  the  elemntary  functions  of  p  and  in  the 
presence  of  the  tabular  values  of  b^a/dr^  is  shown  in  Rabinovich  [l4]. 

It  proved  to  be  insufficient  to  limit  ourselves  to  one  term  in 
the  curvilinear  porticm  when  passing  into  the  temper atwe  region 
below  Boyle ts  temperature.  Therefore  on  seven  isotherms  (T  «  33*26; 
43.16}  53.16}  63.16}  73.l6r  83.16;  95.16)  were  constructed  the  dif¬ 
ferences 

ot* 

Analysis  of  these  differences  on  the  listed  isotherms  showed  that 
in  the  interval  of  change  of  the  cited  densities  o>  *  0  Vo  1  (d^  ■ 

«  0  to  345)  they  may  all  be  determined  by  adding  one  more  meadder  to 
the  cwvillnear  portion  of  the  equation  of  state,  l.e.. 


d*a 

a . 


r? . 


(5) 


In  order  acceptably  to  satisfy  the  initial  values  of  d^c/dr*  in 
the  relatively  slnple  expression  of  the  curvilinear  portion  for  teim- 
peratures  T  >  T^g  the  interval  of  change  in  the  cited  densities  was 
limited  to  <0  «  1,  since  the  substantial  c<mgplications  caused  in  the 
analytical  ejqpression  in  order  to  extend  the  interval  with  respect 
to  <0  are  unjustified. 

WS  set  up  equstions  for  the  two  isotherms  +50  and  <^0^0 
reflecting  the  experimental  data  with  an  accuracy  of  up  to  O.ljt  in 
order  to  determine  the  linear  portion  of  equatim  of  state  (2) . 
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After  the  eubatltutlon  of  the  derived  ejqpreealons  of  the  elemen¬ 
tary  functlona  In  Eq.  (2)  the  equation  of  atate  for  hydrogen  aaaumea 
the  final  form 


«  -  0,48374  «i  +  0.95074  «*  —  0.90968  «.» +  0,91 770  a.*  —  0,19778  + 

+  (1  +  0,23510  »  + 0,01394  «»*  + 0,07029  »*- 0,05824  «>«  + 0,01273  o)‘)i:+ 

+  (— 1,90330  0) —0,77848  «»*  +  0,71 1 10  «*  -  0,66808  +  0,09700  a>*  + 

+  0,02043  ••)  exp  ( —  0, 127  x)  +  ( —  0,65366  «  —  0, 14624  •*  +  3,27170  «>*  — 

-  3,44870  0.*  +  1,18672  0,12090  «)•)  exp(-  1.905x). 

Equation  of  atate  ($)«  true  nthen  T  >  In  the  range  of  the  cited 
denaltlea  of  a>  -  0  to  2.4  (to  ■  3000  kg/cm*)  and  idien  T  <  Tg  In 
the  change  Interval  of  a>  from  0  to  1,  waa  iiaed  to  calculate  the  thermal 
and  caloric  valuea  with  the  aim  of  comparing  them  with  the  ejq;>erlmental 
data. 

Qood  correapondence  of  the  ejqperlmental  and  calculated  valuea  la 
obaerved  In  the  whole  temperatwe  rangej  the  maximum  devlatlona  In 
the  Individual  point  a  do  not  exceed  the  err  or  a  of  the  ejq;>erlfflent  Itaelf. 
The  laotherma  200  and  300°C«  where  the  experimental  valuea  of  o’... 
lie  along  both  aldea  of  the  calculation  ourvea  and  reach  +  0,5»  ere 
exceptlona.  The  dlacrepancy  of  the  experimental  data  m  theae 
laotherma  waa  noted  even  In  analyala  of  the  Initial  valuea. 

On  the  -150^0  and  O^C  laotherma  (Table  1)  the  coincidence  of 
calculated  valuea  of  Internal  energy  with  the  data  In  Nlchela  et  al. 

[13]  la  coagpletely  aatlafactory;  a  aomewhat  poorer  coincidence  of 
valuea  la  obaerved  on  the  130^0  laotherm.  (hily  that  the  130^0  lao- 
therm  waa  the  laat  one  on  the  laoehorea  and  that  In  differentiation 
In  that  work  [13]  Inaccurate  derivative  valxiea  might  have  been  obtained 
can  explain  the  foregoing. 
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Table  2  ccHsparee  the  calculated  values  of  Isochorlc  specific 
heats  with  the  data  In  Michels  et  al.  [15]«  The  average  deviations 
do  not  exceed  2^  and  this  Is  coiapletely  permissible  In  specific  heat 
values . 

The  results  of  the  comparison  enable  us  to  draw  the  following 
conclusion.  Separating  the  curvilinear  portion  of  the  equation  of 
state  by  twofold  differentiation  of  the  Isochores  and  on  this  basis 
determining  the  volumetric  and  temper atxire  fiinctlons  enable  us«  with¬ 
out  deviating  from  the  basic  reqxilrement  that  the  thermal  prcqpoertles 
of  the  substances  ultimately  be  accurately  described,  to  analyze  and 
describe  the  derivatives  which  considerably  Influence  the  caloric  valuee 
and  also  reliably  to  extrapolate  the  equation  of  state  Into  the  high 
temperature  region. 

In  the  calculation  of  the  thermodynamic  values  the  following 
values  of  the  constants  were  adopted! 


Pj^  -  15.225  kg/cm*;  Tj^  -  55.26®  Kj  v^  -  52.248  l/kgj  M  -  2.156; 
420.654  kg-BV^  •  deg;  A  -  426.94"^  kcal/kg-m. , 

Extrapolating  the  Equation  of 
State  Into  the  High-Temperature 
Reglwi 

In  the  hlgh-temperature  reglcxi 
the  thermodynamic  prcqpertles  of 
substances  even  at  high  pressures 


are  chiefly  described  by  the 

Fig.  1.  Isotherms  of  sec<»Ml 

derivatives  of  0*  *  (d^c/^r*)  *  second  and  more  rarely  by  the 
■  f (d.)  for  the  following  tem¬ 
peratures  t  1)  -259  Cj  2)  -220|  third  vlrlal  coefficients,  there- 
5)  -200j  4)  -15OJ  5)  -lOOj  6)  . 

-^1  7)  Oj  8)  lOOj  9)  200;  10)  the  trustworthiness  of  their 

400OC. 

extrapolatl(«i  determines  the 


reliability  of  the  extrapolation  of  the  eqxuitlon  of  state  as  a  whole. 

If  we  collect  all  the  members  of  equation  of  state  (6) which  stand  1 
front  of  <0  we  may  derive  the  analytical  expression  for  the  second 
vlrlal  coefficient 


B  -  0,48374  +  0.235 10 1  - 1 .90337  exp  (-  0. 1 27 1)  -  .  . 

—  0,65566  exp  (—1.905  t).  ' 

The  values  of  the  second  vlrlal  coefficient  In  the  temperature 
range  from  to  -(4000°C  were  calculated  by  this  formula.  In 

order  to  test  the  truth  of  the  values  of  B  In  the  experimental  range 
of  temperature  change  the  data  of  various  auth(»>8  on  the  second  vlrlal 
coefficient  were  plotted  on  the  calculated  curve  B  <■  f  (t)  (Fig.  2) . 

As  Is  evident  from  this  figure  the  calculated  curve  passes  In  the 
very  neatest  way  through  the  points  taken  from  Michels  et  al.  [3]^ 
Holborn  and  Otto  Nljhoff  and  KSesom  [9]*  Michels  and  Qoudeket 
[16].  We  determined  the  values  of  B  at  the  temperatures  of  300  and  400°C 
from  data  on  hydrogen  compressibility  on  those  Isotherms  In  Bartlett 
[1]  and  Wlebe  and  Qaddy  [2]. 

The  laws  of  statistical  physics  were  adduced  to  estimate  the 

* 

reliability  of  the  calculated  values  of  the  second  vlrlal  coefficient 
In  the  temperature  reglm  beyond  the  limits  of  the  experiment. 

Statistical  physics  gives  this  expression  for  the  second  vlrlal 
coefficient: 


2KArf(I— «  *7)r«dr. 
0 


(8) 


If  we  use  the  Lennard-Jones  potential  [6,  12] 
the  reciprocity  between  the  molecules  u(r)  then  we 


as  a  function  of 
may  Integrate 
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TABUS  1 


Conqparlson  of  Calculated  Valws  of  Internal  Energy  AUp 


(iccal/mole)  with  Bata  in  Michele  et  al.  [13]* 


u> 

fm 

,  •# 

— 

ISO 

1  0 

I  ISU 

1 

A  u 

A  u 

A  Up 

1  A  u 

1 

0 

O.OOOCOO  j 

0 

0 

0 

0 

0 

0 

40 

0,115974 

-  9,69 

—  9,60 

—  6,68 

-  7,13 

—  3,68 

-  4.67 

80 

0,231948  i 

—  19,21 

—  19,13 

—  13,07 

—  13,90 

—  6,95 

—  8.83 

120 

0.347922  j 

28,58 

—  28,56 

—  19,22 

-20.33 

-  9.88 

-12,47 

160 

0,463896 

—  37,81 

1—  37,87 

-25,15 

-26,42 

-12,49 

—  15,61 

200 

0.579870  ! 

—  46,88 

—  47.04 

—30,84 

-32, 16 

—14.77 

—  18,2,7 

240 

0,695844  1 

—  55,80 

—  56,04 

1  -36.29 

-37,56 

-16,72 

-20,35 

280 

0.811818 

-  64,55 

1  —  64,86 

-41,45 

1  -42.62 

!  —18,27 

-21.91 

320 

0,927792 

—  73,09 

-  73,48 

—46,28 

-47,32 

-19,38 

-22,87 

360 

1,043766  , 

—  81,37 

-  81,87 

-50,73 

-51;58 

1  -19.95 

-23,17 

400 

1 . 159740 

-  89,26 

-  90,01 

—54,74 

-55,38 

1  —19,92 

-22.74 

440 

1,275714  1 

96,94 

-  97,87 

-58,24 

-58,65 

—  19,20 

-21,56 

480 

1,391688 

-104,27 

-1(K^,41 

-61,18 

-61,33 

-17.71 

—19,56 

520 

1,507662 

-111,18 

-112,59 

—63,48 

-63,38 

-15,34 

-16,71 

560 

1,623636 

-117,56 

-119,32 

—64,99 

-64.71 

—11.95 

—  12.97 

600 

1,739610 

-123,51 

-125,56 

—65,87 

-65,25 

—  7,69 

—  8.33 

640 

1.855584 

—65,82 

-64,84 

—  2,25 

—  2,80 

680 

1,971558 

-64,85 

—63,48 

r  4.34 

4-  3,76 

720 

2,087532 

-62,91 

-61.13 

-rl2,14 

-1-11,43 

760 

2.203506 

-59,89 

-57.70 

-i-21.19 

-i- 20,30 

800 

2,319480 

—55,38 

-53,14 

1  32.02  i 

-j-30,43 

840 

2.435454 

-50,43 

-47,39 

>^43, 11 

-41,89 

-aiiai^wii  . 

1 

1  :0.66 

.43  1 

-rl 

1,98 

Fig.  2.  Bependence  of  second 
▼Irlal  coefficient  B  on  tea^ 
perature  tC^C)  according  tot 
l)  Michels  et  al.  [3]!  2) 
Holborn  and  Otto  [4jj  3)  Nlohels 
et  al.  [16]|  4}  Bartlett  [1]  and 
Vlabe  and  Gaddy  [2jj  3)  Nljhoff 
and  Keeson  [9]. 


and  the  theoretical  values  of 
second  rlrlal  coefficients  and  of 
function  F(x) .  Solid  line)  cal- 
oulatod/  dMhed  line)  theea*etloal. 
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Eq. (8)  using  the  known  H-functlon  and  obtain  a  convenient  forimila  for 
calculating  the  second  vlrlal  coefficient 

=theor  ' 


Where 


6,  =  -^  r.S'rl, 


kT 

u 


(9) 


Tho  function  P(2)  with  a  c<M>pllcated  analytical  form  was  tabulated 
by  Hlrschfelder  [17]  f<a*  different  values  of  argument  £. 

TABIE  2 

Co«q?arlson  of  Calculated  Values  of  Specific  Heats  c^^ 


(kcal/mole  •  deg)  with  Data  In  Michels  et  al.  [13]* 


d 

w 

— 

150 

0 

150 

% 

Cv 

s 

Cv 

Cv 

0 

0.000000 

3,73S 

3,736 

4,852 

4,852 

4,994 

4,994 

40 

0.115074 

3,755 

3,751 

4,873 

4.869 

5,012 

5,013 

SO 

0,231948 

3.773 

3,769 

4,893 

4.886 

5,031 

5,032 

120 

0.347922 

3.793 

3,788 

4,017 

4,904 

5,051 

5,063 

160 

0,463696 

3,812 

3,807 

4,940 

4,920 

5.071 

5,079 

200 

0,579870 

3,831 

3,829 

4,964 

4.941 

5,092 

5,104 

240 

0.60S644 

3,852 

3,850 

4,988 

4,962 

5,113 

5,131 

2S0 

0,611818 

3.872 

3,873 

5,013 

4.985 

5,135 

5,158 

320 

0,927792 

3,886 

3.895 

5,039 

5,009 

5,158 

5,189 

360 

1,043766 

3,913 

3,918 

5,066 

5,033 

5,181 

5,220 

400 

1,150740 

3,930 

3,946 

5,094 

5,059 

5.206 

5,254 

440 

1,275714 

3,953 

3,976 

5,124 

5,086 

5,231 

5,291 

4S0 

1.391688 

3,983 

4,004 

5,155 

5,113 

5,258 

5.330 

520 

1,507662 

4,004 

4,036 

5,187 

5,144 

5.286 

5,369 

560 

1,623636 

4.031 

4,071 

5,221  j 

5,176 

5,316 

5,410 

600 

1,739610 

4,060 

4,112 

5,257  1 

5,210 

5,348 

5.454 

610 

1,855584 

5,294 

5.244 

5,380 

5,494 

660 

1,971558 

5,331 

5,272  I 

5,415 

5,534 

720 

2,087532 

5,874 

5,302 

5.450 

5,578 

760 

2.903506 

5.416 

5,338 

5.487 

5,624 

600 

2,319480 

5,460 

5,374 

5,527 

5,665 

640 

2,435454 

5,503 

5,406 

5,663 

5,711 

anHoi«N«^g 

i0,4 

+0.5 

+  1.2 

Ite  succeeded  In  obtaining  alsmt  cosvlete  coincidence  of  the 
curves  coaipared  not  by  only  parallel  transfer  of  the  coordinate 
system  but  also  by  turning  It,  cp.  Kasavchlnskly  et  al.  [19]«  The 
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angle  of  turning  for  this  case  proved  to  be  very  small  and  did  not 
substantially  ameliorate  the  coincidence  of  the  curves.  In  the  very 
low-ten5)erature  region  of  T  -  20  to  45°  K  the  curves  deviate  con¬ 
siderably.  This  can  be  explained  by  the  fact  that  the  theoretical 
assumptions  get  experimentally  confirmed  for  hardly  any  substances 
at  tenqperatxzpes  close  to  the  critical. 

The  good  agreement  of  the  theoretical  and  calculated  curves  in 
the  experimental  teinperature  range  T  ■  45°  to  675.16°  K  and  also  in 
the  temperature  beyond  the  bounds  of  the  experiment  allows  us  to 
figure  that  the  calculated  values  of  the  second  vlrlal  coefficient  are 
reliable  and  may  be  used  to  deteamlne  the  thermodynamic  properties  of 
hydrogen  at  high  tenqper atures . 

In  conclusion  let  us  note  that  in  the  cited  equation  of  state 
the  detailed  tables  of  the  thermodynamic  properties  of  hydrogen  in 
the  temperature  interval  t  ■  -253°  to  +4000°C  at  in'essures  p  ■  1  to 
500  l(g/cm4  were  calculated  without  allowing  for  dissociation  of 
molecules  at  high  temperature;  and  the  1  —  s  diagram  was  thus  con¬ 
structed. 


Nomenclature 

a^(M),  a,  (<»),Hw)>r  (*»).  ••  functions  of  the  cited  density  coj  ^(t)  Is 

the  mcmotcwilcally  decreasing  functim  of  the  cited  temperature;  N  is 
the  Avogadro  number;  k  is  Boltzmann<s  constant;  u  »  u(r)  is  the 
potential  energy  of  reciprocity  between  two  molecules  at  distance  £ 
frcBi  each  other. 
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Nr.  Copies 


HSADQ,Ui\RTBRS  USAP 
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1 

1 


APSC 
SCPTR 
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25 

5 

1 
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